Pressure dependence of diffusion in simple glasses and supercooled liquids 
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Using molecular dynamics simulation, we have calculated the pressure dependence of the diffusion 
constant in a binary Lennard- Jones Glass. We observe four temperature regimes. The apparent 
activation volume drops from high values in the hot liquid to a plateau value. Near the critical 
temperature of the mode coupling theory it rises steeply, but in the glassy state we find again 
small values, similar to the ones in the liquid. The peak of the activation volume at the critical 
temperature is in agreement with the prediction of mode coupling theory. 
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Diffusion in glasses and their melts has been stud- 
ied intensively for many years. These efforts are stim- 
ulated both by the technological importance of glassy 
and amorphous materials and by the desire to under- 
stand the physics of disordered systems in general and 
the liquid to glass transition in particular. Despite this 
effort there is still no agreement on the nature of diffusion 
on an atomic level or on its change at temperatures near 
the glass transition. This even holds for simple densely 
packed glasses, such as binary metallic ones, see for re- 
cent reviews Refs. |, |, || | 

In a hot liquid, diffusion is by flow, whereas, in the 
glass well below the transition temperature, it will be 
mediated by hopping processes. One key question is the 
transition between the two regimes. For fragile glasses, 
such as most polymers and amorphous metallic glasses, 
mode coupling theory (MCT) predicts an arrest of the 
homogeneous viscous flow in the undercooled melt at a 
temperature Tc, well above the glass transition tempera- 
ture Tg . ^ Hopping processes will suppress the predicted 
singularities and will become the dominant diffusion pro- 
cess near Tc. 

The nature of the hopping process is another issue of 
controversy. Is it by a vacancy mechanism, similar to 
diffusion in the crystalline state, or is it via a collective 
process inherent to the disordered structure ? Investiga- 
tions are hampered by the fact that glasses are thermo- 
dynamically not in equilibrium, and one observes aging. 
The diffusion constant of a glass which has been relaxed 
for a long time will be considerably lower than the one 
of an "as quenched" glass. 

In crystalline materials the pressure dependence of the 
diffusion constant can often be used to identify the dif- 
fusion mechanism. For thermally activated diffusion the 
diffusion constant can be described by an Arrhenius law 



D{T) = Do exp (-ff/fcT) 



(1) 



where Dq is a pre-exponential factor and H is the activa- 
tion enthalpy. In a more component system the diffusion 
constants and the parameters describing them will be 
different for each component. For simplicity's sake we 
drop the indices indicating the component here and in 



the following, whenever possible. 

Using V = dG/dp with G ^ H - TS one obtains the 
activation volume for a diffusion by a single jump process 
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In crystals one finds that the second term is only a mi- 
nor correction which can be neglected. For diffusion me- 
diated by defects, the activation volume splits into two 
terms, T4,ct = + V"™, a formation volume and a 
migration volume V"^ 

For diffusion via thermal vacancies the formation vol- 
ume dominates and V^ct varies between 0.6f2 and Ifi, 
where is the average atomic volume. For the migra- 
tion part one estimates V"^ ~ 0.1 SI. Concomitantly in 
crystals high values of V^ct are taken as a signature of a 
thermally activated diffusion mechanism. 

Assuming that also in the glass the first term in Eq. |^ 
dominates, one usually describes, also in amorphous ma- 
terials, the pressure dependence of diffusion by an appar- 
ent activation volume 
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Experiments on a number of metallic glasses give a 
large spread of values in the range of 0.05 to 1 Low 
values were, e. g., observed in CogiZrig ||^ where no sig- 
nificant isotope effect is observed [||. This result can be 
interpreted in terms of a collective diffusion mechanism 
inherent to the glassy structure. The situation is not so 
clear for the case of large activation volumes T4ct- The 
activation volume of 0.9 SI for diffusion of Zr in Co92Zr8 
was interpreted as indication of diffusion via vacancy like 
defects rather than by an inherent mechanism as in 
the case Co. Values of around 0.5S7 have been observed in 
several materials ^ ^ . Such values are also found 
in materials where the vanishing isotope points to dif- 
fusion by collective jumps. JT^ Whether collectivity can 
induce migrational activation volumes of the order 0.5 to 
1 S7 is still open. 
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In the liquid state the diffusion constant can be fitted 
by a Vogel-Fulcher-Tammann law (VFT) 
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From this we get the activation volume 
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From MCT one derives B alternative expressions 
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In both expressions the diffusion constant extrapolates to 
zero, however, with one crucial difference. In the VFT- 
expression, Eq. |[ this happens at Tq^'^, far below Tg, 
whereas in MCT diffusion vanishes at, T^'~^'^, well above 
Tg. Due to the crossover to the glass the difference be- 
tween the two expressions is normally not sufficient to 
distinguish between them with certainty. This is differ- 
ent for the pressure dependence where both expressions 
give singularities, below and above Tg, respectively. 

There is only one measurement of 14,ct in the under- 
cooled melt. Knorr et aZ.jlJ] find for Ni diffusion in the 
bulk metallic glass melt Zr46.75Ti8.25Cu7.5NiioBe27.5 val- 
ues between 0.35 and 0.64 ft. 

Whereas experiment so far only provides circumstan- 
tial evidence for the diffusion mechanism, molecular dy- 
namics simulation can give both the quantities seen in 
experiment and show the underlying atomic processes. In 
a previous simulation on NiZr the activation volume was 
estimated as 0.36 fl from the change of barrier hight.]!^ 
The underlying process was a collective jump of a chain 
of atoms. 

The aim of the present paper is to present a systematic 
study of the pressure dependence of diffusion as func- 
tion of temperature. In order to relate closely to other 
work the simulations were done for the well studied bi- 
nary Lennard- Jones system (LJ), described by the inter- 
atomic interaction: 

V^JiR) = 4e,, \{a^J/Rf^ ~ {a,j/Rf + Ayi? + B^j 

'(8) 

For the parameters Cij and ct^ we took the commonly 
used values of Kob and Andersen : tAA — ^ ~ <Jaa = 
a — 1, tBB — 0.5, cfBB ~ 0.88, cab = 1-5 and aAB = 0.9. 
Different from these authors we use a larger cutoff radius 
i?c = 3f7. To avoid spurious cutoff effects we introduce, 
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FIG. 1: Diffusion constants (majority ^-atoms, diamonds, 
and minority _B-atoms, circles) and density (triangles) at zero 
pressure against inverse temperature (all in reduced units). 
The dashed and dotted lines show the fits with the MCT and 
VFT expressions, respectively 



similar to the shifted force potential , the parameters 
Aij and Bij to ensure continuity of the potential and 
its first derivative at the cutoff. All masses were set to 
ruj = 1. In the following, we will give all results in the 
reduced units of energy, saa, length, aAA, and atomic 
mass mA- 

We have used this potential previously to calculate the 
diffusional isotope effect parameter in the liquid jl^ and 
the evolution of the dynamic heterogeneity both below 
and above the glass transition From these calcula- 
tions we had well aged samples at different temperatures 
for zero external pressure. 

The calculations were done at each temperature with 
constant volume, corresponding to zero pressure, and pe- 
riodic boundary conditions. The pressure was monitored. 
The time step was At = 0.005. A heat bath was simu- 
lated by comparing the temperature, averaged over 20 
time steps, with the nominal temperature. At each time 
step 1% of the temperature difference was adjusted by 
random additions to the particle velocities. This proce- 
dure assured that existing correlations between the mo- 
tion of atoms were only minimally affected. At each tem- 
perature we had 8 independent samples, each of 5488 
atoms in a ratio 4 : 1 oi A- and i3-atoms. The diffusion 
constants were calculated from the asymptotic slope of 
the atomic mean square displacements 

i?(T) = l lim((R^(t + 0-R'(Ok.' (9) 

where ( )i^t' indicates averaging over all atoms of the 

particular species and over starting times. 

Fig. |l| shows the densities and diffusion constants for 
the zero pressure configurations. From the change in 
slope of the volume expansion we estimate the glass tran- 
sition temperature as Tg « 0.35e//c. The diffusion con- 
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FIG. 2: Calculated activation volume (majority A-atoms, di- 
amonds, and minority B-atoms, circles) versus temperature. 
The open circles represent the experimental data on Ni dif- 
fusion in a bulk glass melt, scaled to LJ-units. The dashed 
and dotted lines are fits with the VFT expression, with one 
or two VFT temperatures, respectively. 



stant can be fitted very well both by the VFT and MCT 
expressions. For both fits we assume that the respective 
transition temperatures are the same for both species. 
For the VFT expression, Eq. |[ we find the parameters 
Tq^ft ^ o.244e/A: for both species and E^^'^ = 1.092 
and 0.986e for A and B atoms, respectively. The val- 
ues for the MCT-expression, Eq. ^, are Tc — 0.36e/fc for 
both species and 7 — 1.87 and 7 = 2.02 for A and B, 
respectively. The two temperatures, Tg and Tc are very 
close to each other, but are much lower than the value 
Tc = 0.435e/A:, reported for simulations at constant den- 
sity p = 1.2. 1^ This reflects the strong dependence of 
the glass transition on density or pressure. We find for 
zero pressure a density of p = 1.15 at T^. 

To calculate the pressure dependence we did additional 
runs with higher and lower densities, respectively. At 
each temperature the zero pressure samples were taken 
as starting configurations which were then compressed 
(expanded) and subsequently aged. The change of den- 
sity was 2% at the highest and 0.5% at the lowest tem- 
perature. The resulting pressure change after aging was 
kept this way within the approximate limits of ±0.2e/(T^. 
Aging varied between 100 000 and several million time 
steps, depending on temperature. 

After calculating the diffusion constants the activation 
volume was evaluated from Eq. ^. Fig. ^ shows the result- 
ing values together with the fit with the VFT-expression, 
Eq. H One can clearly distinguish four temperature re- 
gions. At the highest temperature V^ct ~ 0.6 fi, nearly 
identical for both components. Upon cooling the activa- 
tion volume drops to a plateau value Vact ~ 0.3 fl. The 
larger component has, as one would expect, a slightly 
higher activation volume. Below T = 0.5 the activation 



volume rises sharply for both components and reaches a 
maximum of order f2 around T = 0.4, near Tg and T^^""". 
In the glassy state it drops again to a value below 0.3f2. 

The drop of the activation volume in the liquid corre- 
lates nicely with the drop of the isotope effect parameter 
reported for the same system. [|l8| 

The onset of the increase to the maximum near Tg 
coincides with the onset of a pronounced curvature in 
the Arrhenius plot of D{T), Fig. [I| At this temperature 
and below the isotope effect parameter indicates collec- 
tive motion. |l^ From our previous work and from 
the work on the LJ-system under high pressure ||2^ we 
expect a predominance of chain motion with increasing 
chain length upon cooling. At these temperatures diffu- 
sion exhibits a pronounced dynamic heterogeneity over 
times of several 100 (mCT^/e)^/^.(l9) The time scale of 
this heterogeneity, in turn, correlates with time scale of 
the transition from /3- to a-relaxation. ||l^ 

Our values for the activation energy in the glass are 
of similar magnitude as the estimate of Teichler [|l5| for 
the NiZr system. We find higher activation volumes for 
less well aged samples. In other simulations of metallic 
glasses collective jumps of chains of atoms |2^ 
have been observed, an indication of an inherent collec- 
tive diffusion mechanism. 

The experimental values for Ni diffusion in the bulk 
glass Zr46.75Ti8.25Cu7.5NiioBe27.5 [0 show an increase 
of the activation volume below T = 704 K. When we 
plot (open circles in Fig. ^ the experimental values in LJ 
units by scaling with the temperature where the break in 
the diffusion constant occurs, i. e. without a free pa- 
rameter, they closely follow the LJ result. In particular, 
they show the strong increase of Kct upon cooling to- 
wards T^^"^ . This apparent agreement poses, however, 
some questions. The two temperatures Tg and Tc nearly 
coincide in the LJ system, under the present quench 
and aging times. There is experimental evidence that in 
Zr46.75Ti8.25Cu7.5NiioBe27.5 the two temperatures differ 
strongly: Tg = 592 K and Tc w 875 K.||| According 
to these numbers the increase of the activation volume 
would occur at Tg and not at Tc and, therefore, would 
not be the effect predicted by MCT. The experimental 
evidence is not sufficient to clear this point. There is a 
large experimental uncertainty on the measured values of 
the activation volume, the apparent increase might be an 
artifact of the fitting. Also the MCT critical temperature 
might be overestimated. Clearly better experimental val- 
ues are needed. 

As we have seen above, both VFT and MCT expres- 
sions predict the onset of a divergence in Vact{T). The 
dotted line in Fig. |^ shows a fit with the VFT ex- 
pression, Eq. ||. We assume that the VFT tempera- 
tures are the same for both components, independent 
of pressure. The fit works relatively well. However, 
the resulting parameters are somewhat odd. We find: 
QT^-PT/dp = 0.3 and dE^'^'^ /dp = for both compo- 
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observed in the simulation. This singular behavior can, 
e. g. be used to test whether there is one or several 
critical temperatures in a multicomponent system. 



0.4 0.8 
T ( units of e/k ) 



1.2 



FIG. 3: Calculated activation volume (majority A-atoms, di- 
amonds, and minority B-atoms, circles) versus temperature. 
The dashed line is a fit with the MOT expression using a 
common critical temperature for both components. 



nents and d\nD^^'^/dp = -0.329 and -0.2561 for the 
the A and B-components, respectively. That means, ac- 
cording to the fit there is very little change of the "ac- 
tivation energy" and a strong change of the pre- factor. 
Allowing for different values of dT^^'^ /dp the fit is some- 
what improved. Its main shortcoming, the too slow rise, 
remains. 

The alternative fit with the MCT-expression, Eq. |^, 
is displayed in Fig. ^. As above, we assume a common 
value of T^'^^"^ for both components and all pressures. 
We omitted in the fit the logarithmic term in Eq. |^, 
i.e. we neglected the weak dependence on a possible 
pressure variation of 7. We get an excellent fit (dashed 
line) with dT^'^'^/dp = 0.028 for both components and 
din D^'^^ /dp = -0.338 and -0.278 for A and B atoms, 
respectively. The square deviation for this fit is an order 
of magnitude less than the one using the VFT-expression. 
Contrary to the VFT case it reproduces the sharp uprise. 
We take this as an indication that the divergence of the 
activation volume would in deed be at T^'-^'^ and not at 
yVFT^ Of course the glass transition intervenes before 
the divergence is reached. 

Kob and Andersen report for their constant vol- 
ume simulation a pressure of p « 2.5 at their T^^"^ — 
0.435. Extrapolating our value using the fitted pressure 
derivative we find T^'^'^{p = 2.5) « 0.43 in excellent 
agreement with Kob and Andersen's value. 

In conclusion, we find, both in the liquid and the glass, 
activation volumes of around 0.3 atomic volumes. This 
correlates with a high collectivity seen in the isotope ef- 
fect. In the hot liquid, where diffusion is no longer by 
collective motion, the activation volume rises to about 
0.6 atomic volumes. Cooling towards the critical tem- 
perature mode coupling theory predicts a 1/(T — Tc) sin- 
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